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1 .  INTRODUCTION 


In  a  nonequilibrium  beam-plasma  system,  the  time-average  nonlinear  force 
<F a>  on  a  relativistic  test  particle  due  to  the  total  electromagnetic  field 
consists  of  a  part  <F aa>,  acting  on  the  actual  charge  of  species  a,  and  a  part 
<Fdp >'  acting  on  the  dynamic  polarization  charge  surrounding  the  particle. 
Thus, 


<Fa>  =  cFaa  >  +  <Fdp>  .  (1) 


The  polarization  charge  arises  from  the  redistribution  of  other  particles 
due  to  the  induced  fields  produced  in  the  neighborhood  of  the  particle  as  it 
moves  through  the  beam-plasma.  For  a  relativistic  test  particle  the  time- 
average  force  on  its  actual  charge  (bare  charge)  under  conditions  of  the  Born 
approximation  for  plasma  is  given  by1'* 


<Foa> 


f(1)  +  ?(2) 

*  n  +  F~ 


where 


F(1) 
t  a 


+  Va.Ek 

k  - 

to  +  i6 


and 


►(2) 


dk-j 


6(w  +  u>i  -  (k  +  ki  ).va) 

( u  -f  i  6 )  ( u1  +ifi) 


*  EkiEklj[i?Aif  (k,,k)  +  ^1  Ai  jl)*  (ki  ,k  )  ] 


(2) 


(3) 


(4) 


and 


•a.  V.  Akopyan  and  V.  H.  Tsytovich,  Bremsst  rah  lung  in  a  Nonequilibrium  Plasma,  Fix.  Plazmy,  ± 
(1975),  673  (Sov.  J.  Plasma  Phys.,  ±(197 5),  371). 

*H.  E.  Brandt,  Nonlinear  Force  on  an  Unpolarized  Relativistic  Test  Particle  to  Second-Order  in 
the  Total  Field  in  a  Nonequilibrium  Beam-Plasma  System,  Harry  Diamond  Laboratories,  HDL-PRL-82-7  (May 
1982),  to  be  published  as  HDL-TR-199S. 


Here  t  is  the  time,  ea  is  the  actual  charge  of  the  particle  of  species  a,  ma 

is  its  mass,  k  =  (jc,  u>)  is  a  wave  four-vector,  6  is  a  small  imaginary  part  of 

the  frequency,  <5(x )  is  the  Dirac  delta  function,  E.  .  is  the  Fourier  transform 
til 

of  the  i  component  of  the  total  electric  field,  va  is  the  particle  velocity, 
YQ  =  (l  -  v2/c2)-l/2,  and  c  is  the  speed  of  light.  The  condition  that  equa¬ 
tion  (2)  hold,  namely,  the  Born  approximation,  is  that 


ea|  Ea 


“pe  |  Pa 


where  pa  is  the  particle  momentum  and  uj  is  the  electron  plasma  frequency. 

^  In  the  present  work,  the  time-average  of  the  dynamic  polarization  force 
<Fdp>  is  derived  to  fourth  order  in  the  total  field  for  a  slowly  varying 
nearly  spatially  independent  background  with  no  external  fields.  It  is  shown 
to  be  dependent  on  the  linear,  second-  and  third-order  nonlinear  conductivity 
tensors.  The  result  of  this  calculation  agrees  with  that  of  Akopyan  and 
Tsytovich. 1  It  is  important  in  calculations  of  collective  radiation  processes 
and  the  conditions  for  the  occurrence  of  radiative  instability  in  nonequilib¬ 
rium  beam-plasma  systems. 


2.  THE  DYNAMIC  POLARIZATION  CHARGE  DENSITY 

The  dynamic  polarization  charge  which  surrounds  the  test  particle  results 
from  the  background  distribution  function  being  disturbed  in  the  neighborhood 
of  the  particle.  Thus,  the  associated  current  density  is  given  by 


t  v  fa3Ps  +  f  (S) 

=>dp  =  l  es  /  - -  vs[f 

v  s  J  (2tt)3 


where  the  sum  is  over  all  species  and  f  and  f  '  are  the  perturbed  and 
regular  background  distribution  functions,  respectively,  for  species  s.  It  is 
assumed  that  there  are  no  external  fields  and  the  perturbation  in  the  distri¬ 
bution  is  due  to  the  field  associated  with  the  interaction  between  the  test 


V.  Akopyan  and  V.  N .  Tsytovich,  Breosstrahlung  in  a  Nonequilibrium  Plasma,  Fit,  Plazny,  _i 
(1975),  673  [Sov.  J.  Plasma  Phys.,  1  (1975),  371). 


particle  and  the  other  particles  in  the  system.  The  dynamic  polarization 
charge  density  is  related  to  the  current  in  equation  (7)  by  the  equation  of 
continuity,  namely. 


3pdp 
3t  + 


0 


(8) 


Equation  (8)  assumes  local  conservation  of  the  dynamic  polarization  current. 
The  Fourier  decompositions  of  the  polarization  charge  and  current  densities 
are  given  by 


-  /dk  pdE,k.i <*•*-*>  (91 

and 

3dp  =  /  ifdpkei(**^"I1,t>  »  (10) 

where 


dk  =  d^ic  da) 


(11  ) 


Using  equations  (9)  and  (10)  in  equation  (8),  then 

♦  -»• 

k* jdpk 

p ,  ,  =  - —  . 

dPk  u)  +  i  6 


(12) 


Taking  the  Fourier  transform  of  equation  (7)  one  has  also 


■^dpk 


(s) 


,R(s) 

-  fk  Jv_ 


(13) 


Substituting  equation  (13)  in  equation  (12)  one  obtains  the  following  expres¬ 
sion  for  the  Fourier  transform  of  the  dynamic  polarization  charge  density: 


Pdpk 


k'vs 

a:  +  i5 


(14) 


The  particle  distribution  functions  f^s*  are  determined  by  the  relativistic 
Vlasov  equation,  namely, 


where  the  relativistic  relation  between  velocity  and  momentum  is  given  by 


-► 


and  the  force  Fs 


is  given  by 


+  v, 


VSXBJ 


(16) 


(17) 


Here  E  and  B  are  the  total  electric  and  magnetic  fields.  The  distribution 
function  and  the  fields  can  be  expressed  in  terms  of  their  Fourier  transforms. 
Thus,  for  example. 


f(s)  =  /  dk  f^'e1^^^)  . 


(18) 


In  terms  of  the  Fourier  transforms,  equation  (15)  becomes 

fkS>  =  “7 - — - 77  /  <*1  <*k26(k  -  k!  -  k2)Fsk  -V  f^  .  (19) 

i(w  -  k*v  +  i 6 )  1  ps 

v  s  ‘ 

Expressing  the  distribution  functions  as  power  series  in  the  total  field,  one 
has 


43)  «  f£(s)  +  I  4s)(n) 

n=1 


(20) 


where  f|Ms ^  describes  the  background  for  the  nonequilibrium  beam-plasma. 
The  fk^  describe  the  perturbation  in  the  background  due  to  the  total  elec¬ 


tromagnetic  field.  Assuming  a  slowly  varying,  nearly  spatially  independent 
background  which  in  zeroth  approximation  is  space  and  time  independent  and 
denoted  by  f 0  ^ ,  then 


f£<s>  =  (270-*  Jd3?dt  fR(0)e-i(5t.?-uJt) 


(21  ) 


Performing  the  integration  in  equation  (21),  then 


where  6(k)  =  63(kJ6(u>)  is  the  four-dimensional  Dirac  delta  function.  For 
stationary  turbulence,  equation  (22)  is  exact.  It  follows  from  equations  (19) 
and  (20)  by  iteration  that 
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SK1  ps 
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1 

/  dkidk2  6(k  -  k,  -  k2)Fskl*7p  fk 
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•*  +  1  .*.  +  R(s ) 

x  F  ,  *7  -  F  ,  *7  f, 

sk,  p  .  c  ♦  +  .  £\  sk_  p  k, 

3  s  -  k  *v  +  xo)  5  rs  6 


Substituting  equation  (22)  in  equations  (23,24,25)  and  using  the  properties  of 
the  delta  function  to  do  some  of  the  integrations,  then 


:(s){1) 


i  (u  -  k»vg  +  iS) 


fv4  fR(C 

sk  p  p 
s  s 


.  (s ) ( 2) 


rr - —  /  *,«2  s(*  -  x,  -  k2)?s  .5 
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fR(0) 


-  F  -7  f 
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ft, 


V 


/  dk1dk2d)t3dk4  6(k 


•  8  i^u>2  "  ^2*Vs  +  1°J  3  s  l ( -  kA*vo  +  16 J  4  s 


4  4  s 


The  Fourier  transform  Fsk  of  the  Lorentz  force  Fs  is  given  by 


Fsk  -  eS(4  +  ^s x®k )  • 


However,  from  Maxwell's  equation, 


=  -VxE  , 


it  follows  that 


-»■  1 

=  oj"+  i6  kXE)c  ' 


and  substituting  equation  (31)  in  equation  (29),  then 


Fsk  *  es 


*  SSX(jcxEk) 

*  +  ”  w  +  i6  J  * 


Then  using  the  vector  identity,  one  has 


VgX(j?XEk)  =  (vs«EkK  -  (Vg*£)&k  , 


and  equation  (32)  becomes 


Fsk  =  es 


l*\  ”  «  +  i6/  +  u+ifij 


Next  substituting  equation  (20)  in  equation  (14)  one  obtains  for  the 
Fourier  transform  of  the  dynamic  polarization  charge  density 


l  fd3Ps  S-v, 

n*1  J  (2ir)3  u  +  : 


rs  (s)(n) 
iS  * 


Similarly,  the  Fourier  transform  of  the  dynamic  polarization  current  density 
equation  (13)  becomes 


dpk 


- 1 


i  / 

n=1  J 


d3ps  >  (s ) (n) 
v  f. 


(2ir) 


3  s  k 


Equation  (35)  may  be  written  as 


r  (s) 
Pdpk  -  1  Pdpk  » 


where 


(s)  _  r  is; 
Pdpk  n|1  Pdpk 


(s)(n) 


and 


#  w  ,  r  k-v 

(s)  (n)  =  I  S  _ s 

dpk  J  (2n)3  s  io  +  i 


_  f(s><n) 

6  k 


Similarly,  equation  (36)  may  be  written 


(36) 


(37) 


(38) 


(39) 


■>  r  -»-(s) 

3 dpk  ~  L  3 dpk 
s 


(40) 


where 


+(s) 


3  dpk 


(41  ) 


and 


t(s) (n) 
3  dpk 


♦  (s ) (n) 
esvsfk 


(42) 


Equations  (37)  and  (40)  express  the  Fourier  transform  of  the  dynamic  polar¬ 
ization  charge  density  and  current  density  in  terms  of  a  sum  over  the  contri¬ 
butions  of  each  species.  Equations  (38)  and  (41)  represent  the  latter  as 
expansions  in  the  total  electric  field,  the  nth  order  terms  of  which  are  given 
by  equations  (39)  and  (42).  The  first  three  orders  can  be  obtained  by  substi¬ 
tuting  equations  (26)  to  (28)  and  (34)  in  equations  (39)  and  (42). 


3.  THE  DYNAMIC  POLARIZATION  CURRENT  DENSITY 


The  currents  ?clp^2'/  and  ?ip^3^  are  the  linear,  second-order 

nonlinear,  and  third-order  nonlinear  dynamic  polarization  current  densities, 


respectively.  The  linear  and  nonlinear  electric  field  dependence  is  given  by 
equations  (42),  (26)  to  (28),  and  (34). 


Proceeding  to  reduce  the  linear  current  one  has,  using  equations  (42)  and 
(26), 


+ 

/d  3p  e  v 

_ 1 _ *  * 

.  _  .  3  ,  +  +  ,  Fsk  7] 
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(43) 


Substituting  equation  (34)  in  equation  (43)  then,  equation  (43)  becomes 
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>  R(0) 
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Using  the  Einstein  convention  with  implicit  summation  over  repeated  indices, 
equation  (44)  may  be  rewritten  as 
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Equivalently,  then,  the 
by 


linear  dynamic  polarization  current  density  is  given 
Ddpki1’  =  "ifoOEk*  ,  (46) 


where  the  linear  conductivity  tensor 


is  given  by 
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°i l 
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2V.  N.  Tsytovich,  Theory  of  Turbulent  Plasma,  Consultants  Bureau,  Plenum  Publishing  Corp,  (New 
York,  1977). 


Proceeding  to  reduce  the  second-order  dynamic  polarization  current  using 
equations  (42)  and  (27),  then 
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Substituting  equation  (34)  in  equation  (48),  then 
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Equivalently,  then,  the  second-order  nonlinear  dynamic  polarization  current 
density,  equation  (49),  may  be  rewritten  as 
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where  the  second-order  nonlinear  conductivity  tensor  j  jj(k,ki  ,k2  )  for  species 
s  is  given  by1'3-0 
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Equations  (50)  and  (51)  are  in  complete  accord  with  Tsytovich. 9  He  has  ab¬ 
sorbed  a  factor  of  (-es/o)10)2)  in  equation  (50)  into  his  definition  of  the 
second-order  conductivity.  It  is  to  be  noted,  however,  that  the  Fourier 
transform  and  normalization  conventions  used  by  Tsytovich8  are  identical  to 
those  used  here,  whereas  those  used  by  Akopyan  and  Tsytovich1  are  not.  Note 
also  that  the  tensor  defined  by  their  equations1  (18)  and  (20)  differs  from 
equation  (51)  above  in  that  the  first  complex  denominator  a) 
equation  (51)  here  is  implicitly  a t  -  Jc*vs  -  i6  there.6'7 


K  *vs  +  i  6  in 


Recently,  some  new  exact  symmetries  of  the  second-order  nonlinear  conduc¬ 
tivity  tensor  equation  (51)  have  been  discovered3-8  and  related  to  long- 
established  approximate  symmetries  related  to  the  Manley-Rowe  relations, 
crossing  symmetry,  and  the  nondissipative  nature  of  the  nonlinear  current. 
Also,  a  useful  new  polynomial  representation  for  the  tensor  was  obtained  in 
which  all  derivatives  are  removed  and  the  pole  structure  is  clearly  exhibit¬ 


ed.  3-8  The  symmetries  are  useful  in  the  reduction  of  collective  radiation 


probabilities. 1  Specifically  the  new  exact  symmetries  are  given  by6  8 
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The  approximate  symmetries  follow  from  equations  (52)  and  (53)  when  resonant 
wave-particle  interactions  are  negligible.  For  example,  under  this  condition 
the  well-known  approximate  symmetry — equation  (2.83)  of  Tsytovich2 — may  be 
obtained  from  either  equation  (52)  or  equation  (53).3 * * 6**8'8 

Proceeding  to  reduce  the  third-order  polarization  current  using  equations 
(42)  and  (28)  and  integrating  over  one  of  the  delta  functions,  then 
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Substituting  equation  (34)  in  equation  (54)  and  renaming  wave  vector  variables 
of  integration,  then 
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Equivalently  then,  the  third-order  nonlinear  dynamic  polarization  current 
density — equation  (55) — may  be  written  as 
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Equations  (56)  and  (57)  are  also  in  complete  accord  with  Tsytovich. 9  He  has 
absorbed  a  factor  of  (-eg/ai.,  uijOJj )  into  the  third-order  nonlinear  conductivity. 
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4.  THE  NONLINEAR  DYNAMIC  POLARIZATION  FORCE 

The  dynamic  polarization  force  F^p  on  a  test  particle  is  given  by  the 
Lorentz  force  acting  on  the  polarization  charge  accompanying  the  particle. 
Thus, 

Fdp  =  l  j  <*3r[pd5>  (r,t)!(?,t)  +  j(^p)(?,t)xB(?,t)]  .  (58) 

The  fields  E(r,  t)  and  B[r,  t)  are  the  electric  and  magnetic  fields  induced  by 
the  test  particle  interaction  with  the  other  particles  in  the  system.  It  is 
assumed  that  there  are  no  external  fields.  In  terms  of  the  Fourier  decomposi¬ 
tions,  equation  (58)  becomes  the  following: 

Fdp  =  l  f  d3?dk1dk2  ei^1,?-u,1t)ei^2’r-a)2t) 

r  (s )  -»  >(s )  +  -i  (59) 

x  L  Pdpk  i  Ek  2  +  j  dpk  i  xBk  2  ]  » 

or  performing  the  integral  over  space,  then 
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The  time-average  dynamic  polarization  force  is  given  by 

>  1  ft/2  + 

<FdP>  =  &£  t  J_t/2  FdP(f)df  .  (61) 

Substituting  equation  (60)  in  equation  (61)  and  replacing  the  time  limits  of 
integration  by  the  infinite  limit,  then 
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Performing  the  time  integration,  then  equation  (62)  becomes 
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Using  the  property  of  the  delta  function  to  integrate  over  k1 ,  and  renaming 
the  remaining  integration  variable  k2  to  be  k,  then  equation  (63)  becomes 
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Using  equations  (12)  and  (31),  one  has 
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Next  using  the  vector  identity 
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together  with  equation  (12)  and  the  fact  that  un5(tu)  =  0,  then  equation  (66) 
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Therefore,  substituting  equation  (68)  in  equation  (65)  one  has 
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Substituting  equation  (41)  in  equation  (69),  then 
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Next,  substituting  equation  (46)  in  equation  (71),  then 
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Replacing  the  variable  of  integration  k  by  -k  in  equation  (72),  and  realizing 
that  because  the  fields  are  real 
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then  equation  (72)  becomes 
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An  alternative  form  is  obtained  by  introducing  an  integrated  delta  function  in 
equation  (72)  to  obtain  an  equivalent  expression,  namely, 
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Next  substituting  equation  (50)  in  equation  (71),  and  noting  that  the 
delta  function  is  even,  one  obtains 
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Equation  (76)  is  in  complete  agreement  with  equation  (18)  of  Akopyan  and 
Tsytovich1  since  from  equation  (51)  as  noted  earlier  it  follows  that 


sij  i(“k,ki ,k2  )  -  “Si j fc(k,ki ,k2 )  ,  (77) 

where  j  ^(k,^  ,k2 )  designates  the  conductivity  tensor  appearing  in  Akopyan 
and  Tsytovich*  in  which  evidently  the  first  complex  denominator  is  implicitly 
u  -  5t»vg  -  i 6.  This  may  be  seen  from  equation  (18)  there  where,  because  of 
the  delta  function,  one  has  effectively  “  k2,fc^»k2)»  Assuming  that 


V.  Akopyan  and  V.  N.  Tsytovich,  Bremsst  rah  lung  in  a  Nonequilibrium  Plasma,  Fix,  Plazmy,  1_ 
(197 5),  673  ISOV.  J.  Plasm  Phys.,  l_(1975),  371}. 


equation  (50)  holds  there  also,  except  for  the  differing  Fourier  transform 
convention,  then  because  of  the  delta  function  only  +  k2»k.j,k2)  can 

enter.  Equation  (18)  of  Akopyan  and  Tsytovich1  has  an  additional  factor  of 
(2n)-9(2ir)-3  which  apparently  arises  from  different  Fourier  transform  and 
background  normalization  conventions.  For  example,  in  their  equation  (5)1  the 
Fourier  transform  convention  employed  has  a  factor  of  (2 it)  ^  in  the  inverse 
Fourier  transform  in  the  integration  over  the  three-dimensional  wave  vector 
space,  and  a  factor  of  1  for  the  integration  over  frequency,  giving  a  total 
factor  of  (2ir)~3,  whereas  here  a  total  factor  of  1  is  used  as  in  equation  (9), 
for  example.  Also,  the  Fourier  transform  itself  has  a  factor  of  (2it)  1  there 
and  (2TT)"4  here  as  in  equation  (21),  for  example.  There  is  also  another  addi¬ 
tional  factor  of  ( 2 n ) ~ 3  in  their  equation  (18)  1  which  is  apparently  due  to 
differing  background  normalization.  Because  of  the  different  Fourier  transform 
convention  alone,  the  counterpart  of  their  equation  (18)  1  would  have  an  addi¬ 
tional  factor  of  ( 2 rr ) 3.  Apparently  it  has  been  absorbed  into  the  normal¬ 
ization  of  fp|0)  there.  In  short,  the  fp|0)  there  must  be  (2 n)  3  times  that 
here.  Alternatively  if  the  normalization  is  in  fact  the  same  as  that  here, 
then  there  is  an  erroneous  factor  of  (2ir)  ^  appearing  there.  Also  the  factor 
of  1/6  in  Akopyan  and  Tsytovich1  arises  from  the  explicit  symmetrization 
chosen  there. 

Next  substituting  equation  (56)  in  equation  (71)  and  noting  that  the  delta 
function  is  even,  one  obtains 
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It  is  to  be  noted  that  equations  (78)  and  (57)  are  in  apparent  agreement  with 
equations  (19)  and  (21)  of  Akopyan  and  Tsytovich.1  Evidently  in  the  conduc¬ 
tivity  tensor  T?  j  jjn,  (k,k.,  'k2'k3  )  9iven  by  ^e^uation  (21)  there,1  the  first  ^two 
complex  denominators  are  implicitly  uj  -  jt»vg  -  i  <$  and  uj  +  +  1?1  )*vg  - 

i6,  respectively,  and  also  there  is  no  overall  factor  of  i  as  there  is  in 
equation  (57)  here.  Therefore 


ZijL(~k'k1'k2'k3) 


“iTijL(k'k1  ,k2'k3 ) 


lA.  V.  Akopyan  and  V.  I*.  Tsytovich,  Bremsstrahlung  in  a  Nonequilibrium  Plasma,  Fiz.  Plazmy,  _i 
(1975),  673  (Sov.  J.  Plasma  Pbys.,  1  (1975),  371). 


Thus,  comparing  equations  (78)  and  (79)  with  equation  (19)  of  Akopyan  and 
Tsytovich, 1  one  finds  that  they  are  in  complete  agreement.  Again  the 
additional  factors  of  (2u)_12  and  (2n)“3  and  1/24  in  their  equation  (19)  1 
apparently  arise  from  differing  Fourier  transform  and  background  normalization 
conventions  and  explicit  symmetrization,  respectively. 

In  summary,  then,  using  equation  (70),  the  time-average  dynamic  polar¬ 
ization  force  on  a  relativistic  test  particle  to  fourth  order,  in  the  field  in 
a  nonequilibrium  beam-plasma  system  for  a  slowly  varying,  nearly  spatially  in¬ 
dependent  background  with  no  external  fields,  is  given  by 


•*■(0) 

<Fdp>  =  Fdp 


?(1)  +  ?(2> 

Fdp  +  Fdp 


(80) 


The  linear,  second-order  nonlinear,  and  third-order  nonlinear  dynamic  polar- 
ization  forces  F^p',  F^p',  and  are  given  by  equations  (75),  (76),  and 

(78),  respectively. 


5.  CONCLUSION 


»  - 


An  expression — equations  (80),  (75),  (76),  and  (78) — has  been  obtained  for 
the  time-averaged  dynamic  polarization  force  on  a  relativistic  test  particle 
to  fourth  order  in  the  total  field  in  a  nonequilibrium  beam-plasma  system  for 
a  slowly  varying,  nearly  spatially  independent  background  with  no  external 
fields.  This  result  has  been  used  in  the  work  of  Akopyan  and  Tsytovich  in  the 
theory  of  collective  bremsstrahlung  in  nonequilibrium  plasmas. 


The  present  work  together  with  related  work  by  the  author 10' *' + ' *  is 
important  for  ongoing  work  in  calculating  collective  radiation  processes  and 
conditions  for  the  occurrence  of  radiative  instability  in  relativistic  beam- 
plasma  systems. 
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